Thermodynamics of the Double-Layer Quantum Hall Systems 
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In this paper we apply the exact solution of the sine-Gordon model to describe thermodynamic 
properties of the soliton liquid in the incommensurate phase of the double-layer quantum Hall 
systems. In this way we include thermal fluctuations and extend to finite temperatures the results 
obtained by C.B. Hanna, A.H. MacDonald and S.M. Girvin [Phys. Rev. B 63, 125305 (2001)]. In 
addition we calculate the specific heat of the system. While the results obtained for the sine-Gordon 
model are available in a temperature interval {Q,Tc), where Tc = 8nps, ps the pseudospin stiffness, 
they can be applied in the bilayer system up to temperatures STbkt, where Tbkt = 7rPs/2 is the 
vortex mediated Berezinskii-Kosterlitz-Thouless transition temperature. Above this temperature 
the operators cos l3ip and cos(27n9//3) are both relevant and the system is in a phase with coexisting 
order parameters, is the dual field of and /3 is the sine-Gordon coupling constant. We provide 
numerical estimates for thermodynamic quantities for the range of parameters relevant for GaAs. 

PACS No: 71.10.Pm, 65.40.-f, 65.50.-Fm 



Strong electron-electron interactions in Quantum Hall 
Double Layer systems lead to interesting collective 
effects 0. Perhaps, the most remarkable among them 
is the interlayer coherence which is established when the 
interlayer spacing d is of the order of the magnetic length. 
In that case the Coulomb interaction between electrons 
on different layers is comparable to the interaction on the 
same layer [s!] . The dynamics of this coherent state is con- 
veniently described in terms of the pseudospin variable 
for which the up and down states refer to the electron 
being on the upper or the lower layer whereas real 
spins are totally frozen. 

At finite interlayer distance the symmetry in pseu- 
dospin space is lowered down to U(l) such that the cor- 
responding Ginzburg-Landau free energy is of an easy- 
plane anisotropic ferromagnet. The latter symmetry is 
further broken by interlayer tunneling. In the limit of 
strong anisotropy one can treat pseudospins as belonging 
to a plane such that at low temperatures the free energy 
is simplified down to the classical sine-Gordon model [3] : 
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where (0 — Qr) is the phase angle describing pseu- 
dospins in the XY plane at different positions. Q = 
(27r(i/0o) B|| X z is the parallel magnetic field wave vec- 
tor, with (/)o — hc/e being one flux quantum between the 
layers, t is the tunneling energy which generally depends 
on magnetic field Q as well as on 



t = toe 



(2) 



niz = vi — V2 is the layer imbalance and ps is the pseu- 
dospin stiffness which depends on as well 



Ps arises from the loss of the Coulomb exchange in the 
presence of a phase gradient (in what follows we neglect 
the dependence of ps on m^, considering that the layer 
imbalance is small or equivalently vi = 1 — V2 ~ 1/2)- 
In the above formulas / is the magnetic length equal to 
{hc/eB±y/^ , where B±^ is the strength of the magnetic 
field perpendicular to the layers. 

At higher temperatures Eq. ^ does not properly de- 
scribe the pseudospin system as it does not take into ac- 
count the 2tt periodicity of the field (p. This free energy 
should be supplemented by the contribution originating 
from the vortex configurations of angle variable (p: 



d xcos2n9, — = , (4) 
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where 5*0 is the thermodynamic action of a vortex and 
is the dual field of (j>. The scaling dimensions of this term 
and the tunneling term are respectively 
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such that 



d^dt = 1/4 
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In what follows we shall apply non-perturbative results 
derived for the sine-Gordon model to the bilayer systems. 
Strictly speaking, this approach is valid only below the 
Berezinskii-Kosterlitz-Thouless transition temperature 



Tbkt = 



TPs 
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Ps = (l- ml)pE 
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when dv > 2. We shall first discuss the region T < 
Tbkt and then briefly discuss what happens at higher 
temperatures. 

To study two-dimensional classical sine-Gordon (SG) 
model ([T]) we use the well known analogy with the quan- 
tum (l+l)-dimcnsional SG model, at T = 0. The action 
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of the quantum SG model corresponding to ([T]) is: 



where 
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where /S^ = T/ps, ip = <j)/f3, = t/^AirPT). The field 
|h| = 27r|Q|//3^, coupled to the sohton topological charge 
density, plays the role of a chemical potential. The two 
models have the same partition function. The current 
field theory description is valid only at distances much 
greater than the magnetic length I which serves as the 
ultraviolet cutoff. 

The ground state energy Eq of the quantum model is 
related to the free energy of the classical model: 
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This is a fundamental formula which we use to establish 
a link between the exact solution of the quantum sine- 
Gordon model and the model describing Quantum Hall 
double layer. 

The SG model, given by Eq. ([T]), contains two com- 
peting periodicities: the periodic potential tends to lock 
the system into a commensurate configuration with the 
field (j) being locked to a minima of the cosine potential, 
whereas the gradient term prefers the field configuration 
he (j) = Qx. This competition takes place only for tem- 
peratures smaller than the critical temperature 
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below which the cosine potential is relevant. Notice that 
Tc is 16 times larger than Tbkt and therefore the 
interlayer tunneling is highly relevant in the entire area 
of validity of the sine-Gordon description. In this area 
there is a critical value of chemical potential Qc{T) such 
that for Q > Qc, the competition is resolved by a forma- 
tion of a liquid of solitons (domain walls). Each soliton 
interpolates between minima of the cosine potential. The 
transition into such incommensurate state occurs when 
the soliton energy equals the chemical potential. 

According to [1| , numerical values of the parameters for 
a typical GaAs double layer QH sample are the follow- 
ing. The effective mass is m* « 0.07me, total particle 
aerial density is tit — 1-0 x 10^^ cm^^, layer (midwell 
to midwell) separation is d = 20 nm and tunneling en- 
ergy is to = 0.1 meV. For such a sample we would have 
I « 12.6 nm, d/l « 1.6, huic ~ 6.9 meV for i^t = 1, and 
e^/iirel w 8.8 meV. This gives ps » 0.03 meV, which 
corresponds to Tbkt ~ 0.5K. 

At Q = Qc a single soliton is introduced in the system. 
The value of Qc as determined by the exact solution ^] 
is given by the equation 
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In the limit T — > this equation coincides with the ex- 
pression used by Hanna et al. [1] : 
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Qc{T) is a monotonously decreasing function of T , as 
shown on Fig. 1. The maximal value of Qc is achieved 
for T = 0; substituting the numerical values used for 
GaAs, we obtain the estimate for a critical magnetic field 
necessary to observe the soliton liquid: 

Qc(T = 0)/«0.3 , B,(T = 0) -O.l^o/W . (14) 

From this estimate we also see that since the value of 
{Qcl/'2) < 0.15 is always small, it is not always necessary 
to take into account the field dependence of the tunneling 
integral in Eq. (|11[) . (This term is however important in 
cases of strong magnetic fields and mistakenly it is not 
accounted in pj] in all asymptotic cases Q — *■ -|-cxd they 
considered. The correction terms to quantities like M|[, 
Xsoi etc., resulting from the Q dependence of t, should 
decrease exponentially when Q — > +oo instead of the 
power low decay ofjlj]). 

The authors of [l|] used T ^ limit to describe the 
soliton state. To establish whether such description can 
be extended to finite T, we have to recall some funda- 
mental facts about the sine-Gordon model. 

The particle spectrum of the quantum SG model con- 
sists of solitons and antisolitons of mass Ms and for 
T < 1/2 also of their bound states 
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Recall that the sine-Gordon description is valid below 
Tbkt which corresponds to t < 1/16. In the case the 
system will be exposed to an external "magnetic field" 
only the solitons' and antisolitons' energy will be affected. 
The breathers have zero charge and do not interact with 
external magnetic fields. The solitons acquire additional 
energy —h whereas the antisolitons h, and in the ground 
state only solitons can appear. The soliton's excitation 
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FIG. 1: Plot of 4r A/sol on temperature (r — T/Tc) and mag- 
netic field Q dependence is shown, to is the tunneling energy 
and Tc is the critical temperature above which the soliton 
generating operator becomes irrelevant. The line on the sur- 
face shows the critical magnetic field, above which the system 
crosses to the incommensurate phase. (On the vertical axis 
we show tMsoi instead of the diverging soliton mass Msoi). 
As can be seen, for fixed r, the mass of the solitons decreases 
exponentially as Q increases. The Q and Ms axes have units 
l/l, (the inverse magnetic length, I = 12.6 nm). 



spectrum (in the presence of the field h coupled to the 
soliton charge), is given by the following equations 



EiO) = e{9) = [eo{0) - h] 



+ B 



d0'G{e - e')e{e') (is) 



p{e) = 27r / Ae'a{e') ^Po{e) 



+ B 



<:\e'd{e - e')a{e') . 

-B 

e{6) is the dressed energy whereas eo(^) = Ms coahO 
and Po{0) — Ms sinhfl are the bare energy and momen- 
tum. iy{9) is the rapidity distribution function (in ra- 
pidity space). In the above equation the boundary B is 
defined by the condition e(±_B) = 0, whereas the kernel 
G{9 — 9') is the derivative with 9 of the soliton-soliton 
scattering matrix Sss{9 — 9'). The Fourier transform of 
G{9-9') is given in Rcf. [5]. 

e{9) will have a negative part only ii h > he — Ms, 
where Ms is the soliton mass. This constitutes also the 
condition set on h in order the system to acquire a single 
soliton. Therefore the relationship between the soliton 
mass and the critical field is: 



Qc^^rMsiQc) 



(19) 



In the parameter range relevant for GaAs we have /i(0) ~ 
0.03. 

In Fig. [1] we plot tMs as function of r and Q. The 
soliton mass Ms diverges as 1/r as r ^ 0. On the surface 
of the three dimensional plot (t,Q,4tMs) we show also 
the critical field Qc{t). Its dependence on temperature is 
shown also in the inset of Fig.[T] The vertical dotted line 



in it separates the temperature intervals < T < Tbkt 
and Tbkt <T <Tc. 

At Q 7^ only the solitons' and antisolitons' energies 
are affected. The breathers have zero topological charge 
and do not interact with external magnetic fields. 

As we have said, the sine-Gordon model description is 
rigorously valid only below Tbkt! however one can use it 
as a good approximation even above Tbkt provided the 
correlation length generated by the interlayer tunneling is 
shorter than the correlation length generated by vortices: 



Mi(T) » [exp(-5o)]i/(2^ 



(20) 



where Mi is the mass of the first breather determined 
by Eq. fTSl) . Assuming that 5o ~ 2Tips/T which is the 
thermodynamic action of a skyrmion in the 0(3) non- 
linear sigma model, we estimate that this inequality is 
reasonably fulfilled up to T « 3Tbkt- 

There are two experimentally measurable quantities 
one can extract from the sine-Gordon thermodynamics: 
the specific heat and the magnetization. Both quantities 
are related to the free energy; at t < 1/2 the latter one 
contains two contributions: 
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F2 being related to the ground state energy of the soli- 
tonic Fermi sea. (The solitons have a relativistic disper- 
sion relation and usually their spectrum is parametrized 
in terms of the rapidity 0, in the form pq = M^sinh^, 
eo(^) = Ms coshf?, for the bare values of energy and mo- 
menta. B and the renormalized value of e(9) can be found 
in [11). Therefore Tj 7^ only at Q > Qc- In the com- 
mensurate phase there will be also a contribution from 
the the first part of Ti due to the Q-dependence of Mg'. 



^ j«soi + (QV2)-^cot 

4 1 — T 



2(1 -r) 
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{Mslf (2,4) 
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with Usoi being the soliton density. In the limiting case 
T ^ the subcritical contribution to the magnetization 
is given by 

^M„(<3<ft;T.0).f , m 

where Qc is given by Eq.([T 
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The magnetization 





FIG. 2: Plot of the magnetization on both the C and I phases 
of the double-layer QH system as a function of the magnetic 
field Q and temperature (r = T /Tc) is shown. In the C phase 
the magnetization increases with Q until it reaches the C-I 
transition critical point. The presence of the soliton conden- 
sate in the I phase results in a fast decrease of the magneti- 
zation. 



FIG. 3: The total specific heat of the double layer QH system 
for different values of the magnetic field Q. The presence 
of the soliton condensate leads to a decrease in the specific 
heat, except close to the critical line where it diverges. As Q 
increases the critical point goes towards T = point and the 
divergence changes form from l/\/e to l/[e ln(l/e)]. Here the 
tunneling amplitude is t = 0.1 meV. 



On Fig. 2 we give the dependence of the magnetization 
on the magnetic field and temperature. As can be also 
observed in the figure the presence of the solitons leads 
to a decrease of the magnetization of the system. It takes 
biggest value at small temperatures and close to the criti- 
cal line. For high fields the magnetization vanishes. Close 
to the critical line, in the I phase the thermal fluctuations 
of the solitons renormalize the magnetization decrease by 
changing its behavior from l/ln(l/e) as T — > to ^/e at 
finite temperatures, where e — [t/tc{Q) — 1]. 

In the following we give also the asymptotic value of 
the density of solitons (or the soliton lattice wave number 
Qs = 2ix jLs = 27rnsoi, is the distance of solitons from 
each other) . This is related to the soliton lattice magneti- 
zation by M|| = (47r^ps/(/)o)nsoh where ngoi is the density 
of the solitons in the system in the IC phase and (/)o is 
the magnetic flux quantum. The density of solitons can 
be found as: 



"sol 
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'271 dQ 



' dh 
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h — 2-kQ/(x' is the field coupled to the soliton charge, or 
the number of particles, and — ^/ Ps: [see Eq. ([1])]. At 
finite temperature we find 



[Q , T ^ 0] 



.1/2 



QQc 



16 
3\/2 



1/2, 



8\/2 



3 1 



G(0)e 



(28) 



The last term is a result of the fact that the tunneling 
amplitude is field dependent. 



Similarly along the fixed Q line, Qs, to first order in 
(ctY^^, has square root behavior in T 



Qs = 




(29) 



For smaller temperatures (fixed r) the square root de- 
pendence holds only in a decreasingly small interval Q 
above Qc- At these temperatures for larger Q the depen- 
dence of 1/Ls on Q changes Using the expression for the 
soliton free energy Q we find 



[Q -> Qc , (T -> 0)] 



0. 
Q 



QQc 
Q,QP 



(30) 



ln(l/e) 



where Ms would be the soliton mass in the case the tun- 
neling amplitude is independent of the magnetic field Q. 

On Fig. 3 we give the plot of the specific heat for ^ = 
0.03. At small T in the C phase the specific heat per unit 
area is linear in temperature: 



Cv/T = r^iJ-zi InV 
8 



(31) 



For jji = 0.03 we have 7?^ « 0.14i?. 

In the I phase the soliton condensate leads to a decrease 
of the specific heat, except close to the critical line where 
it diverges. It is worth mentioning that the leading order 
term on the condensate contribution to the specific heat 
changes form from l/[eln(l/e)] when T — > to l/^/e at 
finite temperatures. 
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FIG. 4: Plot of the temperature and field dependence of the 
compressional stiffness tensor element Ku is shown. Close to 
Tc, Kii approaches the constant value pa. As T ^ Tkt, Ku 
reaches quickly values close to the asymptotic value, ps, and 
does not vary much with the magnetic field. 



In the following we calculate the compressional stiff- 
ness element Ku of the stiffness tensor Kij. The plot 
which is given below applies in the case of the bilayers 
up to temperature T « 3Tkt, and for other systems ana- 
logues to SG, in the whole temperature region. We calcu- 
late Ku as the change of the free energy by varying the 
spacing between solitons. This is achieved by calculating 
the second derivative with respect to the component Q^i 
for fixed Q 



Ku{Q,t) 



sol 



(32) 



where /soi is the soliton contribution to the free energy 
and is equal to F2. Substituting 



Qs 



1 dU 



Ps dQ 

one finds for Ku the following expression 

1 /solg/solQQQ 



(33) 
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JsolQQ 
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{dQs/dQ) 
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Qs (dH 



{dQs/dQy 



(34) 



Note that in Ref. 1] the second term is not taken into 
account. 

We calculate Ku numerically and give in Fig. 2] its 
dependence as a function of temperature and magnetic 
field. 

In the following we give the asymptotic behavior of 
Ku, keeping only the leading order term. 

At finite temperatures, near the critical line, Ku in- 
creases as a square root function of the magnetic field 



.1/2 



■ 4V2tM2 



(35) 



As it can be seen on Fig. [H the compressional elas- 
tic constant characterizing the soliton lattice, reaches 
quicker its asymptotic value at higher temperatures than 
at lower ones. At temperatures close to Tc, Ku has val- 
ues very close to ps and changes very little for variations 
of the magnetic field. 

For temperatures towards zero the square root holds 
in a decreasingly small interval Q above Qc- For larger 
Q, outside this interval, the dependence of Ku on Q 
changes. At zero temperatures the Q interval of the 
square root behavior vanishes and the Q dependence of 
Ku is 



[Q-^Qc,T~*0] 
Ql 



' 87r2r2M2 



(36) 



The large Q asymptotics at zero temperature is given 



by 



[g ^ +00, 0] 



Ku ^Ps 1+ (37) 



'64 I 2 Q ) e(Q0V2 



The asymptotics for the case of the tunneling amplitude 
independent of the magnetic field is obtained by substi- 
tuting 1 = 0. 

The large Q asymptotics for finite temperatures and 
field independent tunneling is given by 

[Q^+oo,r^O] (38) 
'18- 84r+ 128r2 - Gir^) 



Ku = P. 
where 

a{T) = 



l + a{T)- 



r(r) r(5/2-T) (?fi*)4(i— ) 



(39) 



r(-r) r(l/2 + r) (3 - 2r)(2T - 1) 
It is not the soliton mass itself, but the combination 't'I 



= 2V^- 
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i2(l-r)J 
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[2(l-r)J 
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r(i + T) 
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which enters as a scale in the expansion above. 

The case with the tunnehng dependent on the mag- 
netic field Kii has the following asymptotic behavior 



1 + a(T) 



e('Q)V2 



(40) 



This is consistent with (|37|) . if in (j40| r = is substi- 
tuted. In (|40p only the highest order term is kept. 

In this work we made use of the exact solution of the 
sine-Gordon model to extend the previous results for the 
bilaycr QH system Q to finite temperatures by includ- 
ing the thermal fluctuations of the soliton lattice. We 
calculated the magnetization and its temperauture de- 
pendence in both C and I phases. In addition we calcu- 
lated the specific heat of the bilayer system which could 
not be achieved with the method of the previous work 
[H . The presence of the condensate results in a decrease 
of both of the specific heat of the system (except close to 
the critical line where it diverges) and its magnetization. 

We observe however that this method cannot be ap- 
plied in the whole temperature interval where the re- 
sults for the sine-Gordon model are available. Above the 
Berezinskii-Kosterlitz-Thouless transition temperature, 
in the interval (Tbkt,Tc) where cos(3(p and cos(27ri?//3) 
are both relevant, the system is in a new phase with coex- 
isting order parameters. (As we argued we could extend 
the interval of SG applicability, with decreasing accuracy, 
up to 3Tbkt)- In this phase this system can no longer 
be described by Eq. ([T]). Furthermore, as pointed out by 
[§| , at higher temperatures the role of the spin degree of 
freedom becomes important and the model used would 



need to be modified. 

We notice finally that different authors have found, 
instead of Tc, the Berezinskii-Kosterlitz-Thoouless tran- 
sition temperature as the point where the soliton lattice 
melts down. 
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